Superconductivity in the two-dimensional t-J model 
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Using computational techniques, it is shown that pairing is a robust property of hole doped anti- 
ferromagnetic (AF) insulators. In one dimension (ID) and for two-leg ladder systems, a BCS-like 
variational wave function with long-bond spin-singlets and a Jastrow factor provides an accurate 
representation of the ground state of the t-J model, even though strong quantum fluctuations destroy 
the off-diagonal superconducting (SC) long-range order in this case. However, in two dimensions 
(2D) it is argued - and numerically confirmed using several techniques, especially quantum Monte 
Carlo (QMC) - that quantum fluctuations are not strong enough to suppress superconductivity. 
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The nature of high temperature superconductors re- 
mains an important unsolved problem in condensed mat- 
ter physics. Strong electronic correlations are widely be- 
lieved to be crucial for the understanding of these ma- 
terials. Among the several proposed theories are those 
where antiferromagnetism induces pairing in the d x 2_ y 2 
channel 0] . These approaches include the following two 
classes: (i) theories based on Resonant Valence Bond 
(RVB) wave functions, with electrons paired in long spin 
singlets in all possible arrangements d||, and (ii) theo- 
ries based on two-hole d x 2_ y 2 bound states at infinitesi- 
mal doping, formed to minimize the damage of individual 
holes to the AF order parameter, which condense at finite 
pair density into a superconductor j|. However, recent 
density matrix renormalization group (DMRG) calcula- 
tions have seriously questioned these approaches since 
non-SC striped ground states were reported for realistic 
couplings and densities in the t-J model Clearly to 
make progress in the understanding of copper oxides, the 
2D t-J model ground state must be fully understood, to 
distinguish among the many proposals. 

In this paper, using a variety of powerful numerical 
techniques, the properties of the t-J model are investi- 
gated. Our main result is that in the realistic regime of 
couplings the 2D t-J model supports a d x 2_ y 2 SC ground 
state, confirming theories of Cu-oxides based on AF cor- 
relations. The t-J model used here is 
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where Cj )Cr =Ci )Cr (1 — n>i^), (...) stands for nearest- 
neighbor sites, and and Si are the electron density 
and spin at site i, respectively. Our study focuses on 
the low hole-doping region of chains, two-leg ladders, 
and square clusters, using different numerical techniques: 
QMC (pure variational and fixed-node (FN) approxi- 
mations), DMRG, and Lanczos. Within our QMC ap- 



proach, it is possible to further improve the variational 
and FN accuracy by applying a few (p < 2) Lanczos 
steps to the variational (p = 0) wave function I'Fy), 
I*p> = C 1 + Efe=i a k H k )\^ v ) §. Non-variational es- 
timates of energy and correlation functions can also be 
extracted with the variance-extrapolation method || . 
Our BCS variational wave function is defined as 
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where Vjq is the projector onto the subspace of N parti- 
cles, Vg is the Gutzwiller projector, which forbids dou- 
bly occupied sites, and J = exp(JV Vijhihj) is a Jas- 
trow factor, defined in terms of the hole density at site 
i, with hi = (1 — — a n( i v i.j being varia- 

tional parameters. The Fourier transform of the pair- 
in g amplitude, f jj, satisfies (|]: f k = A k /(e k — fi + 
y/ (e k — /i) 2 + A 2 ), where e k is the free electron disper- 
sion, \i is the chemical potential, and A k is the BCS 
SC gap function. In particular, it is known that for 
square lattices and periodic boundary conditions (PBC), 
the most relevant contribution to A k has d x 2 —y 2 sym- 
metry, namely A k = A(cosfc 2: — cosk y ) ||. For lad- 
ders, Afc = A x cos k x + A y cos k y was used (the opti- 
mized A x and A y have opposite signs), whereas in the 
ID case A k = Aicos/c + A3COs3fc. In the following, 
VMC denotes results obtained with \^ v ), VMC+pLS 
(with p = 1,2) those obtained with \^ p ), and FN and 
FN+pLS results obtained using the FN approximation 
with |\fv) and \ty p ) as guiding wave function, respec- 
tively. Finally, variance indicates results obtained with 
the variance extrapolation method. 

The wave function Eq. (j^) describes preformed elec- 
tron pairs, expected to become SC within the RVB sce- 
nario ^|J3). An important component of Eq.(|^) is the 
Gutzwiller projector, which at half-filling freezes the 
charge dynamics, establishing quasi-long-range AF or- 
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der |Q. This shows that the projected BCS wave func- 
tion describes magnetic regimes as well. In addition, the 
SC order parameter is not simply related to the pair am- 
plitude A, as in weak-coupling BCS. In fact, at low hole- 
doping, it is proportional to the number of holes, and 
not to the number of electrons. This result is natural 
in hole-pairing theories 0], where superconductivity at 
half-filling is not possible, suggesting that such theories 
maybe similar to the RVB approach if the latter incor- 
porates long-range singlets, and no-double occupancy is 
enforced. Moreover, it was observed that the variational 
parameter A decreases with increasing hole doping ||, 
suggesting a relation of this quantity with the pseudogap 
of underdoped cuprates ||. In hole-pair based theories, 
a similar result is obtained with the hole binding energy, 
finite even at half- filling, playing the role of A M. 
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FIG. 1. S z (q) and energies for (a) ID and (b) two-leg lad- 
ders at zero hole doping. The numerically exact result (open 
circles) was obtained with the FN method. Pairing correla- 
tions for the t-J model in (c) ID (absolute value) and (d) on 
a tilted 2D cluster for parallel (upper points) and orthogonal 
(lower points) singlets. 

To show that the wave function Eq. (^) accurately de- 
scribes not just a SC state but several magnetic sys- 
tems as well, consider the half-filled model on chains 
and ladders. In the first case, the ground state is quasi- 
antifcrromagnetically ordered, with zero staggered mag- 
netization and power-law spin correlations, while in the 
second case there is a finite spin gap in the spectrum, 
and exponentially decreasing spin correlations |l0[. The 



spin structure factor S z (q) = 1/L J2i j e 



iq(Ri-Hj, 



shows a cusp at q = n in ID, and a broad maximum 
at q — (ir, 7r) for two-leg ladders. These features are re- 
markably well reproduced by our variational wave func- 
tion (Fig. 0a-b), which generates robust AF correlations 
at short distances. The state used is also surprisingly 
accurate in energy compared to work based on a pure 
Gutzwiller wave function H, and on numerical exact 



studies. It seems that Gutzwiller projecting the BCS 
wave function allows for a quantitative description of AF 
correlations in low-dimensional systems |ll| . 

Since undoped systems with short-range AF correla- 
tions appear properly described by the wave function 
Eq. (|^), consider now hole-doped AF systems where su- 
perconductivity should emerge according to some the- 
ories Ilj-H]. For this purpose the pairing correlation 



function A^'^(r)- 



Si+ r U S] „ was studied. Here Sj 



(cj -|cJ +M ^ — cj i cJ + „ |) creates an electron singlet pair in 
the neighboring sites + fi). Off-diagonal long-range 
order (ODLRO) is implied if P d = 21im r ^ co y/\A^"(r)\ 
remains finite in the thermodynamic limit. In ID, 
ODLRO is suppressed by quantum fluctuations, but 
A^' v (r) is finite at short distances and the accuracy of 
the wave function Eq.(§) can be assessed, the FN pro- 
viding exact results in ID. As shown in Fig. |l|c, the ID 
pairing correlations are indeed non-zero, although rapidly 
decaying with distance. Our variational wave function 
reproduces accurately the pairing correlations, but only 
when a long-range Jastrow factor is included, otherwise 
the tendency to pairing is overemphasized. The accu- 
racy of the variational wave function is excellent also for 
small 2D clusters, where the exact solution can be ob- 
tained by Lanczos. In this case the quantum fluctuations 
appear not strong enough to destroy superconductivity, 
see Fig. 
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FIG. 2. (a) Average rung hole density and (b) pairing 
correlations for a 30 x 2 ladder, with 12 holes, J—t and OBC. 
In (b) the distance r runs from the center of the ladder. 

Now we consider doped two-leg ladders. In Fig. ^, re- 
sults on a 30 x 2 ladder with open boundary conditions 
(OBC) and J=t are presented [^2|. Our results show 
that in this case the ground state is qualitatively well de- 
scribed by a projected BCS wave function with a density 
Jastrow term, although charge oscillations induced by 
hole-pair formation are strong in DMRG (here considered 
as exact), but considerably weaker with VMC. However, 
it is remarkable that the FN approach gives the correct 
rung density profile, showing that, even in a non-trivial 
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system such as a two-leg ladder with hole pairs, the cor- 
relation functions can be well-controlled by QMC meth- 
ods. In general, it is observed that whenever the VMC 
method is not quantitatively accurate, the proper corre- 
lation functions arc obtained by applying the FN approx- 
imation [^3|. Therefore, this QMC approach combining 
VMC and FN methods represents a novel powerful tool 
to assess the reliability of a variational state , and to 
obtain accurate properties of t-J models for cuprates Q . 
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6x6 cluster FN does not provide qualitative changes in 
the pairing correlations with respect to the VMC out- 
come, and the FN+2LS and VMC+2LS correlations are 
close. In addition, the m=1200 DMRG results lead to 
similar correlations. The agreement among the many 
methods suggests that the QMC method correctly repro- 
duces ground state properties also in 2D systems, where 
the Jastrow term does not suppress the ODLRO present 
in the VMC, providing sizable long-range pairing. 

The most natural boundaries without spurious symme- 
try breaking are PBC, since the finite-cluster eigenstates 
have all the lattice symmetries. Therefore, here the sub- 
sequent effort building toward the main result of the pa- 
per focuses on PBC clusters. For these boundary con- 
ditions, our improved variational calculation is accurate 
and no sign of static stripes has been found at the cou- 
plings investigated, although dynamical effects are still 
possible (see below) . For 8 holes on the 8x8 lattice and 
J/t=0A, our best variational energy per site (FN+2LS) 
is E=— 0.66672(6)t, close to the zero variance extrap- 
olation E=— 0.671(1)2;, and much better than the pure 
variational calculation, E—~ 0.64266(8)i. 



FIG. 3. Pairing correlations (a) and total energy (b) for 
various techniques on 6x6 lattices, m is the number of states 
in the DMRG approach. In (a) the distance r runs along the 
center of the cylinder in the periodic direction. 

After testing the numerical methods, let us now ad- 
dress the main subject of the paper, i.e. the possibility 
of SC in the 2D t-J model. For this case no exact solution 
(analytical or numerical) is available and, therefore, it is 
crucial to perform a careful computational analysis, com- 
paring the results of different techniques. While DMRG 
allows for an almost exact ground-state characterization 
for ID systems and two-leg ladders, unfortunately in 2D 
the results appear to depend on the boundary condi- 
tions . In order to compare the performances of QMC 
and DMRG (with m states kept), not only the standard 
PBC have been considered, for which reliable DMRG cal- 
culations on large 2D clusters presently are not possi- 
ble, but also the cylindrical boundary conditions (CBC), 
open (periodic) in the x (y) direction, and OBC in both 
directions. To reduce the number of variational param- 
eters with CBC and OBC, the Jastrow factor was re- 
stricted to depend only on the distance between sites, i.e., 
v i,j = v i-j- Site dependent chemical potentials /ii were 
also used as additional variational parameters to consider 
possible non-uniform charge distributions. To test the 
stability of our variational wave function directly in 2D, 
a 6x6 lattice with 6 holes was considered. In Fig. [|, a 
comparison among different numerical techniques for sev- 
eral boundary conditions is shown. The large-™ extrap- 
olated energies of DMRG are remarkably similar to those 
of QMC for CBC and OBC, but not for PBC where QMC 
produces substantially better energies. As in ID, for the 




FIG. 4. Superconducting order parameter Pd vs 8 at 
J/t=0.4, using the techniques and cluster sizes shown. In- 
set shows S z (q) at optimal density. 

The main result of our paper is in Fig. |], where the 
SC order parameter Pd is shown for clusters of 64 sites 
at several densities and 242 sites close to optimal dop- 
ing. These results, stable among the many methods used 
here and with weak size effects, are indicative of a ro- 
bust SC ground state away from half-filling in the t-J 
model. In view of the success of the present QMC tech- 
niques to reproduce known results for chains, ladders, 
and OBC 2D clusters, it is reasonable to consider the 
data in Fig. ^ as accurate. In addition, the results are 
in good qualitative agreement with experimental tenden- 
cies, including an optimal doping at (5^0.18. The spin 
structure factor in Fig. ^ (inset) has a broad peak at 
(7r,7r), concomitant with exponentially decaying AF cor- 
relations in real space. A slight tendency toward spin 
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incommensurability (SI) appears in the FN+2LS results 
at optimal density. At other densities 5, such as 0.25, a 
similar mild tendency toward SI was also observed. Then, 
it is conceivable that the SC state may contain very weak 
dynamical stripe tendencies, as recently observed in the 
spin-fermion model p5| . Alternatively, band effects or 
AF correlations across-holes |l(| could be responsible for 
the SI structure. In this doping region, the charge struc- 
ture factor N(q) is found to be basically featureless. 

In addition, a low-density of 8 holes on a tilted 242- 
site cluster was also analyzed to investigate coexistent 
antiferromagnetism and superconductivity, as seen in re- 
cent experiments for underdoped YBCO 0|. The pure 
variational approach shows a small SC order parameter, 
and a vanishing small AF order (Fig. ||). Indeed, at half- 
filling, the projected d x 2_ y 2 BCS state underestimates 
the magnetic order since it has zero magnetization with 
a logarithmically divergent S z (ir, tt) @], and hole doping 
reduces further the AF correlations. Remarkably, the FN 
approach enhances both the SC and AF tendencies. In 
particular, the spin correlations show robust long-range 
order implying that antiferromagnetism survives a small 
density range [jl8],[l9|, as in experiments [jl7|. The FN 
energies (not shown) for 0,8,24 and 40 holes in 242 sites 
are clearly stable against phase separation, but are quite 
close to it [Eoi. 
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FIG. 5. Spin-spin correlations (a) and pairing (b) versus 
distance. In the inset of (b) are the hole-hole correlations. 
The cluster has 242 sites, 8 holes, and J=0.4t. 

Moreover, hole-hole correlations display sharp peaks 
at hole distance (see Fig. ||), indicating the tendency 
to have short-range hole-pairs in the low hole-density 
ground stat e |4| . These pairs were identified in previ- 
ous studies [pi| , illustrating the good agreement among 
different techniques. Using the present QMC method the 
PBC 6x6 two- holes ground state was also analyzed. The 
results (not shown) once again indicated a clear signal at 
distance which always emerges upon Jastrow and/or 
FN improvement from an RVB wave function, where elec- 
trons rather than holes naively appear to be paired. This 
result illustrates the flexibility of the QMC approach, 



and unveils unexpected similarities between the Jastrow 
corrected long-singlet RVB |||| and the hole-pair ap- 
proaches 0, both producing closely related results. 

In conclusion, robust indications of superconductiv- 
ity have been found in the 2D t-J model. The results 
were obtained with an RVB wave function with Jastrow 
factors, flexible enough to reproduce DMRG results for 
chains and ladders and to produce tight hole pairs at low 
hole density, in agreement with other techniques. The 
present effort substantially improves on previous calcu- 
lations where d-wave superconductivity was predicted to 
exist at intermediate J/t H[|, and highlights the power 
of the recently developed QMC methods ||. 
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